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ABSTRACT 


Radon  source  potential  characterization  of  sites  in  terms  of  soil  index 
properties  which  do  not  vary  with  transient  conditions  such  as  moisture  content, 
barometric  pressure,  temperature  and  wind  speed  is  studied.  The  invariant  index 
properties  which  were  found  to  be  critical  for  site  characterization  are  radium 
activity  concentration  in  the  soil,  in-place  dry  density,  porosity,  and  dry  gas 
permeability.  These  properties  can  be  measured  in  situ  or  in  the  laboratory,  or 
estimated  on  the  basis  of  other  soil  index  properties  such  as  grain-size 
distribution  and  Atterberg  limits.  Various  expressions  for  radon  source 
potential  are  reviewed  and  a new  expression  is  formulated  on  the  basis  of  data 
from  areas  of  deep  glacial  terrace  deposits.  Site  exploration  methods  proposed 
include  use  of  the  Standard  Penetration  Test  together  with  a laboratory 
determination  of  radlxim  activity  concentration,  and  a rapid  field  measurement 
procedure  using  a portable  gamma  ray  spectrometer,  a portable  nuclear  moisture- 
density  meter  and  retrieval  of  a soil  sample  for  laboratory  determination  of 
particle-size  distribution.  A plan  to  develop  exploration  protocols,  test  the 
effectiveness  of  the  source  potential  prediction,  and  prepare  a draft 
exploration  standard  is  proposed. 


Keywords:  convective  flow; 
permeability;  porous  media; 
potential . 
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LIST  OF  SYMBOLS 


radium  activity  concentration  in  the  soil  (Bq-kg’^) 

C = radon  concentration  in  soil  gas  (Bq-m'^) 

Cm  = measured  radon  concentration  in  soil  gas  (pCi-l”^) 

^max“  maximum  possible  radon  concentration  in  the  soil  gas 
^60 

Cu  = -j—  = uniformity  coefficient  of  particle  size  distribution 

*^10 

c = coefficient  in  the  Hazen  formula 

D = bulk  diffusivity  coefficient  in  soil  matrix  (m^-s"^) 

Dg  = effective  (interstitial)  diffusion  coefficient  (m^-s”^) 
d = geometric  mean  particle  diameter  (by  size)(m) 

dj^g  = the  particle  size  which  is  equal  or  larger  than  10%  of  the  particles  by 
weight  or  "effective  particle  size"  (m) 
dgg  = the  particle  size  which  is  equal  or  larger  than  60%  of  the  particles  by 
weight  (m) 

e = void  ratio  (dimensionless) 

F = radon  source  potential  (Nazaroff  6e  Sextro,  Bq-s”^) 
f = Emanation  Coefficient  (dimensionless) 

G = Volumetric  Radon  Generation  in  Soil  Pores  (Bq • m” ^ • s ” ^ )* 

H = height  below  soil  surface  (m.) 
h = air  change  period  (h) 

J = radon  flux  (Bq-s'^-m“^) 

k = coefficient  of  gas  permeability  (m^) 

km  = measured  coefficient  of  gas  permeability  in  (m^) 

kw  = Darcy's  water  permeability  coefficient  (m-s’^) 

L = length  of  the  perimeter  crack  (m) 
i = diffusion  length  (m) 

m = degree  of  saturation  as  a fraction  of  pore  space  (dimensionless) 

M = mean  migration  distance  (Tanner,  m) 
n = porosity  of  soil  (dimensionless) 

Pq  = mean  atmospheric  pressure  (Pa) 


*lBq  corresponds  to  the  emission  of  1 radon  atom  per  second  (Ig  of  radium  emits 
3.7-101°  Bq  or  1 Cu) . 
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p = disturbance  pressure  (Pa) 

Pjj  = total  disturbance  pressure  drop  (Pa) 

R = radon  index  (Scott) 

RIN  = radon  index,  (Eaton  and  Scott,  PCi/i) 

RINjj  = radon  index  (Kuntz,  Pci/i) 

r = half  the  width  of  a perimeter  crack  (m  or  mm) 

S = geometric  standard  deviation  of  particle  diameter  distribution 
t = time  (s) 

t3=  time  required  to  displace  fraction  a of  the  total  quantity  of  gas  (min) 

V = seepage  velocity  vector  of  soil  gas  averaged  over  gross  soil 

cross  section  (m-s"^) 

Q = rate  of  flow  (m^-s"^) 

V = radon  index  (Yokel,  dimensionless) 

7^  = dry  density  of  soil  (kg-m‘^) 

7w  = unit  weight  of  water  (Pa-m*^) 

ARn=  decay  constant  of  radon  (2.1- 10'®  s~^) 

fj.  = dynamic  viscosity  of  soil  gas  (Pa-s))  @ 18°  & 760mm  Hg,  /i  = 1,827- 10"^ 
Pa-s,  p(gas)  = 1.204  kg-m”^ 

/iw  = dynamic  viscosity  of  water  (Pa-s) 

Pg  = density  of  soil  grains  (normally  2.65  Mg-m’^) 

V = gradient  operator  (m'^) 

V-  = divergence  operator  (m~^) 

Vp  = the  disturbance  pressure  gradient  (Pa-m”^) 

= Laplace  operator  (m'^) 
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EXECUTIVE  SUMMARY 


There  is  a need  to  identify  and  measure  site  characteristics  which  contribute  to 
the  emission  of  radon  gas  into  the  interior  of  buildings.  In  present  practice 
this  is  usually  attempted  by  insitu  measurements  of  gas  permeability  and  radon 
activity  concentrations  in  the  soil  gas.  These  measurements  are  sensitive  to 
transient  conditions  prevailing  at  any  point  in  time:  the  gas  permeability  is 
sensitive  to  soil  moisture  content;  and  the  radon  activity  concentration  in  the 
soil  gas  is  sensitive  to  variations  in  soil  moisture  content,  atmospheric 
pressure,  temperature  gradients  and  wind  speed.  However,  there  are  measurable 
invariant  soil  index  properties  (independent  of  moisture  content  and  atmospheric 
conditions)  which  are  related  to  the  propensity  of  sites  to  generate  and 
transport  radon  gas  over  time.  This  report  deals  with  the  site  characterization 
for  radon  source  potential  using  these  invariant  soil  index  properties.  The 
advantage  of  this  approach  is  that  the  measurements  performed  in  a site 
exploration  are  repeatable,  a necessary  pre-condition  for  any  measurement 
associated  with  a standard,  and  that  the  calculated  source  potential  represents 
the  propensity  for  radon  generation,  rather  than  the  radon  generation  under  a 
specific  set  of  conditions,  which  may  be  transient.  The  site  exploration 
procedures  considered  utilize  established  soil  exploration  and  classification 
techniques  to  the  maximum  extent  possible. 

The  most  important  invariant  soil  properties  identified  are  radium  activity 
concentration  (Aj^^)  , in-place  dry  density  and  porosity  of  the  subsoil  (7^^,  n)  , 
and  the  dry  gas  permeability  coefficient  of  the  subsoil  (k)  . These  quantities 
either  can  be  measured  or  estimated  on  the  basis  of  soil  index  properties. 
Quantities  , 7^^  and  n can  be  measured  directly  by  various  field  and 
laboratory  techniques.  The  permeability,  k,  can  be  estimated  on  the  basis  of 
grain  size  distribution,  which  in  turn  can  be  measured  by  established 
procedures.  An  insitu  value  of  k (not  corrected  for  moisture  content)  is 
presently  measured  by  air  Injection  or  soil  gas  extraction.  This  measurement  is 
strongly  dependent  on  moisture  content.  In  most  instances  a correction  for  the 
moisture  content  would  have  to  be  made  if  this  measurement  is  used  to  determine 
dry  gas  permeability.  While  equations  which  correlate  soil  moisture  content 
with  gas  permeability  have  been  developed,  presently  available  information  does 
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not  provide  a reliable  basis  for  a moisture-content  correction  of  gas 
permeability  without  the  aid  of  additional,  sophisticated  laboratory  tests. 
Alternately,  the  coefficient  of  gas  permeability  can  be  derived  from  that  for 
saturated  water  permeability,  which  in  turn  can  be  measured  in  the  field  and  in 
the  laboratory. 

Several  expressions  for  radon  source  potential  have  been  proposed  (Equations  28, 
29,  30,  32,  and  33).  An  additional  expression  (equation  31)  is  proposed  in  this 
paper  on  the  basis  of  data  from  New  York  State  (Kuntz  et  al.,  1988).  The 
dimensional  correlation  between  the  permeability  and  radon  entry  depends  on 
whether  the  radon  transport  by  convective  flow  is  dominated  by  transient  or  by 
steady-state  flow  conditions.  Several  of  the  proposed  expressions  are  derived 
from  steady-state  flow  conditions.  However  available  data  indicate  that  the 
radon  concentration  is  approximately  proportional  to  k°-^  (for  instance  Kuntz  et 
al.  1988,  Scott,  1985),  a dimensional  correlation  which  is  closer  to  that 
predicted  by  transient  flow  conditions.  There  are  insufficient  data  to  draw 
conclusions,  because  it  is  difficult  to  separate  the  effect  of  k from  that  of 
other  variables. 

Some  of  the  proposed  expressions  for  radon  source  potential  use  invariant  soil 
properties,  while  others  are  linked  with  site  exploration  procedures  which 
measure  in-place  gas  permeability  and  radon  activity  concentration  in  the  soil 
gas  without  consideration  of  transient  conditions.  All  these  expressions  could 
be  modified  to  utilize  measured  invariant  soil  properties  with  a proper 
allowance  for  the  effect  of  soil  moisture  content. 

Two  exploration  procedures  are  proposed:  (1)  exploration  by  traditional  test 
boring  methods,  such  as  the  Standard  Penetration  Test,  supplemented  by  7- 
spectroscopy  to  determine  Aj^^  and  by  a procedure  to  identify  finegrained  soils 
which  are  likely  to  develop  fissures;  and  (2)  a rapid  field  exploration 
procedure  using  a portable  spectrometer  to  determine  Aj^^ , a portable  nuclear 
moisture-density  gage  to  determine  7^  and  natural  moisture  content,  and 
retrieval  of  a small  soil  sample  for  laboratory  particle-size  analysis  and 
indirect  determination  of  k.  It  is  suggested  that  this  latter  field  procedure 
could  be  combined  with  conventional  percolation  tests  to  provide  more  accurate 
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information  on  permeability.  Whether  method  (1)  or  (2)  is  used  would  depend  on 
the  scale  of  the  exploration  (one  building  or  many  buildings)  and  on  the  need 
for  geotechnical  data  for  other  purposes  (i.e.  method  (1)  would  be  used  if  soil 
borings  are  also  needed  for  other  purposes) . 

A plan  for  further  studies  is  proposed  in  order  to  develop  exploration 
protocols,  test  the  validity  of  radon  source  potential  predictions,  and  prepare 
a draft  standard. 
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1.  INTRODUCTION 


Soil  and  rock  have  been  identified  as  the  most  important  source  of  radon  to 
which  people  are  exposed  [Nazaroff  et  al.,  1988].  It  is,  therefore,  important  to 
identify  the  site  characteristics  which  contribute  to  the  transport  of  radon  gas 
into  the  interior  of  buildings.  The  indoor  radon  level  in  a building  depends  on 
the  potential  of  the  site  to  emit  radon,  and  on  the  characteristics  of  the 
building  and  its  HVAC  system.  Thus  a true  measure  of  the  radon  source  potential 
for  a building  must  consider  site,  as  well  as  building  characteristics.  This 
report  deals  only  with  site  characterization. 

Presently  used  site  characterization  methods  employ  various  techniques  to 
measure  two  quantities:  the  radon  concentration  in  the  soil  gas;  and  the  gas 
permeability  of  the  soil  on  the  building  site  [Tanner,  1988;  Sextro  et  al.; 
1988;  Kunz  et  al.,  1988].  The  measurement  of  these  quantities  is  sensitive  to 
transient  conditions.  The  radon  activity  concentration  in  the  soil  gas  at  any 
point  in  time  depends  on  the  radon  flux,  which  in  turn  is  influenced  by 
variations  and  the  rate  of  change  in  barometric  pressure,  soil  moisture  content, 
wind  velocity,  and  temperature.  In  addition,  the  radon  activity  concentration  in 
the  soil  gas  is  also  sensitive  to  barometric  pressure  changes,  rainfall,  and 
temperature  variations  prior  to  its  measurement  [i.e.  Schumann  et  al.,  1988]. 
The  gas  permeability  of  the  soil  varies  with  the  moisture  content,  which  in  turn 
varies  seasonally  and  could  be  changed  by  future  site  development. 

The  objective  of  this  study  is  to  assess  the  possibility  of  developing  a site 
characterization  procedure  which  is  based  on  measurements  which  are  not 
sensitive  to  transient  conditions.  Consideration  is  also  given  to  the 
development  of  procedures  which  will  utilize  established  soil  exploration  and 
classification  techniques  to  the  maximum  extent  possible.  Equations  for  radon 
production  and  transport  are  reviewed  in  Section  2 and  some  solutions  are 
presented.  The  variables  in  these  equations  and  their  measurement  are  discussed 
in  Section  3.  In  Section  4 various  expressions  for  radon  source  potential  are 
reviewed  and  a new  expression  is  proposed  on  the  basis  of  data  from  New-York 
State.  Options  for  site  exploration  are  discussed  in  Section  5.  A plan  for 
standard  development  is  presented  in  Section  6. 
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2.  RADON  PRODUCTION  AND  TRANSPORT  IN  SOILS 


2.1  Radon  Production  in  Soils 

The  volumetric  radon  generation  rate  in  the  soil  pores  can  be  determined  by  the 
following  equation  [Nazaroff  et  al.,  1988,  p,85]: 

G = ^ ....(1) 

where:  G = volumetric  radon  generation  in  soil  pores  (Bq-m”^s"^)* 

f = emanation  coefficient  (dimensionless) 

Ajj3=  radium  activity  concentration  in  the  soil  (Bq-kgl'^) 

Pg  = density  of  soil  grains  (normally  2.65  Mg-m'^) 
decay  constant  of  radon  (2.1-10'®  s'^) 
n = porosity  of  soil  (dimensionless) . 

*lBq  corresponds  to  the  emission  of  1 radon  atom  per  second  (Ig  of  radium  emits 
3.7-10^°  Bq  or  1 Cu) . 

Equation  (1)  can  be  re-written  as: 

G = (f7dAR.^R„)  ^ ■■■(2) 

where:  7^  = dry  density  of  soil  (kg-m'®) 

In  the  absence  of  any  radon  flux  (say  at  a very  great  depth)  the  rate  of 
volumetric  radon  generation  equals  the  rate  of  radon  decay: 

~ ^Rn*^  ...  (3) 

Therefore : 

C = C,„=  ^ = (f7<,AR.)  ^ ...(4) 

^Rn 

where : 

C = radon  activity  concentration  in  soil  gas  (Bq.m'®) 

^max“  maximum  possible  radon  activity  concentration  in  the  soil  gas 
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2 . 2 Radon  Transport 


The  following  general  equation  can  be  written  for  the  radon  concentration  field 
in  the  soil  (averaged  over  the  total  soil  volume)  [Oliveira  Lourero,  1987]: 

1^  (Cn)  = V-  (DVC)  - V-  (vO)  + n(G  - A^^C)  ...(5) 

where:  V = gradient  operator  (m‘^) 

V-  = divergence  operator  (m'^) 

D = bulk  diffusivity  coefficient  in  soil  matrix  (m^-s”^) 

V = seepage  velocity  vector  of  soil  gas  averaged  over  gross  soil 
cross  section  (m-s'^) 

Eq  (5)  is  for  the  transient  case,  where  there  is  a spaclal  variation  of  radon 
activity  concentration  with  time.  For  steady-state  conditions: 

^=0  ...(6) 

The  first  term  on  the  right  side  of  Eq  (5)  deals  with  diffusive  transport  which 
follows  Pick's  law  and  can  be  modeled  by  the  following  equation: 

= DV^C  - n(V„  - G)  ...(7) 

A solution  for  the  steady-state  diffusion  from  a semi- inf inite  homogeneous  soil 
mass  with  an  activity  concentration  C = ^ great  depth  and  an  activity 

concentration  C = 0 at  the  surface  is  given  by  the  following  equation  for  C at 
depth  z below  the  surface  [Clements,  1974]: 

C(z)  = C,,Jl-exp(-  f)]  ...(8a) 

where:  i = (-^)°-^ 

The  radon  flux  at  the  surface  resulting  from  the  special  distribution  of  C in  Eq 
(8a)  is: 

~ ^^Rn-^^max  • • • (^h) 
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where  H is  known  as  the  "diffusion  length"  and  is  typically  1 m,  A typical  value 
for  Dq  = D/n,  the  "interstitial"  (or  "effective")  diffusion  coefficient  for 
radon  in  soils  of  low  moisture  content  is  on  the  order  of  10”®  m^s'^  [Nazaroff 
and  Nero,  1988].  In  saturated  soil  the  coefficient  may  be  reduced  to  the  much 
lower  value  of  2*10"^°  m^s"^.  An  equation  predicting  the  diffusion  coefficient 
as  a function  of  porosity  (n)  and  saturation  (m)  has  been  proposed  by  Rogers  and 

Nielson,  1988:  Dg  = 7- 10'®  exp[-4(m  - mn^-  m^)]  in  m^  s'^  ...(9) 

The  latter  equation  yields  a value  of  7- 10”®  m^-s”^  for  dry  soils.  In  the  same 
reference  another,  modified  equation  is  proposed  for  highly  compacted  soils.  Eq 
(9)  is  empirical  and  was  obtained  by  fitting  experimental  data. 

The  second  term  on  the  right  side  of  Eq  (5)  deals  with  convective  flow  of  soil 
gas,  the  most  important  transport  mechanism  contributing  to  the  entry  of  radon 
gas  into  buildings.  Within  the  range  of  pressure  gradients  and  soil  particle 
size  distributions  which  is  of  practical  interest  in  conjunction  with  the  radon 
transport  problem,  it  is  reasonable  to  assume  that  convective  soil-gas  transport 
can  be  modeled  by  Darcy's  law  [Darcy,  1856],  which  is  valid  when  the  flow  is 
laminar  and  the  flow  resistance  is  frictional.  The  upper  limit  to  the  validity 
of  Darcy's  law  would  occur  at  high  Reynold's  numbers  where  the  laminar  flow 
breaks  down  and  becomes  turbulent.  This  limit  has  been  observed  to  be  exceeded 
in  gas  flow  through  soils  at  relatively  high  induced  pressure  gradients 
[Matthews  et  al.,  1988].  However,  since  Reynold's  Number  is  directly 
proportional  to  flow  velocity,  and  velocities  are  small  when  the  induced 
pressure  gradients  are  small,  the  laminar  flow  limit  is  not  likely  to  be 
exceeded  in  the  flow  conditions  under  consideration  even  in  gravelly  soils. 
There  is  some  evidence  [Hansbo,  1960]  that  for  the  flow  of  pore  water  there  is  a 
lower  limit  when  particle  sizes  are  very  small  (some  clays)  and  surface 
chemistry,  rather  then  friction  provides  most  of  the  flow  resistance.  This 
range,  if  applicable  to  the  flow  of  soil  gas,  is  also  not  of  interest,  because 
flow  velocities  in  very  finegrained  soils  are  very  small  and  the  dominant 
mechanism  of  radon  transport  is  diffusion. 

For  flow  of  gases  through  a porous  medium  Darcy's  law  can  be  written  as: 

V = -(-)  Vp  . . .(10) 
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where ; 


k = coefficient  of  gas  permeability  (m^) 

H = dynamic  viscosity  of  soil  gas  (Pa’s)) 

(3  18°  & 760mm  Hg,  n = 1. 827-10'^  Pa-s,  p(gas)  = 1.204  kg-m”^ 

Vp  = the  disturbance  pressure  gradient  (Pa-m"^) 

The  concept  of  disturbance  pressure  used  in  Eq  (10)  is  similar  to  the  concept  of 
"excess  pore  water  pressures"  used  in  geotechnical  engineering.  Consequently, 
disturbance  pressure  is  defined  as  the  difference  between  the  pressure 
prevailing  at  any  point  in  the  soil  and  the  atmospheric  pressure  for  the  point 
(which  increases  with  depth) . 

If  it  is  assumed  that  the  disturbance  pressures  in  the  soil  gas  are  small  in 
comparison  with  the  prevailing  atmospheric  pressure,  that  there  are  no 
significant  temperature  variations  within  the  space  considered  which  could 
influence  soil  gas  pressures,  and  that  the  soil  is  homogeneous  and  isotropic 
with  respect  to  permeability,  the  following  general  equation  applies  to 
pressure-gradient  induced  convective  flow  of  soil  gas  [Fukuda,  1955]; 

kPo  3p 

■ V2p  = ~ ...  (11) 

pn  3t 

This  equation  applies  to  the  transient  condition,  where  the  pressure 
distribution  changes  with  time.  For  steady-state  flow  conditions  (when  pressure 
differentials  prevail  long  enough  for  the  gradients  to  stabilize) , this  equation 
becomes : 


V^p  = 0,  or  when  k is  not  constant:  V- (kVp)  = 0 ...(12) 

Two  solutions  to  Eq  (11)  are  examined  below,  Schery  et  al.,  1984,  present  the 
"steady-state"  (stabilized)  solution  to  a sinusoidal  pressure  pulse: 

if:  p = PgSinajp  and  p^  « 

then;  = 0,t)  = Po(^^  )°-^sin(wt  - 135°)  . . . (13) 

where;  = mean  atmospheric  pressure 

p = disturbance  pressure 
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Another  case  of  interest  is  an  atmospheric  disturbance  pressure,  p^ , applied  at 
time  t = 0 at  the  surface  of  a soil  layer  of  infinite  depth.  The  disturbance 
pressures  and  convective  flow  velocities  (or  convective  flux)  as  a function  of 
depth  z at  any  time,  t,  can  be  derived  from  a solution  to  the  general  heat 
equation  presented  by  Carslaw  and  Jaeger,  1959,  p 60,  equations  (5)  to  (8): 

= erf[  I (^)“  =-t-»  5]  ...(14) 

Po  ^ *^^0 

V2(z  = 0,t)  = Vq  = ...(15) 

and  v^(z,t)  = Vg  • exp(-|^|-^)  ...(16) 

When  a layer  of  depth  H resting  on  an  impervious  base  is  substituted  for  the 
infinite  depth  in  the  previous  solution,  the  following  expressions  are  derived 
from  an  approximate  solution  of  the  one -dimensional  consolidation  equation 
[Terzaghi,  1942].  If  it  is  assumed  that  at  constant  temperature  the  volume  of 
pore  gas  is  inversely  proportional  to  the  pressure  (Boyle's  law),  then  the 
volume  of  gas  displaced  by  a disturbance  pressure  p^  from  a column  of  soil  of 
unit  cross  section  area  and  height  H is: 

nHpo 

= “5 — > Po  « Pq 

and  for  the  following  conditions: 

If:  t < 0 ; P = Pq  for  ail  values  of  z 

0 < t < 0.413  and  z = 0 ; P = P^  - p. 

kPq  ° ° 

then:  = 0,r)dr  ==  -2p„(^^^  (17) 

Differentiating:  '^^(z  = 0,  t)  - -pg( — )0-5.^-o.5  ....(15) 

Equation  17  is  for  the  total  radon  flux  from  t=0  to  t.  The  negative  sign  in 
Eqs.(15)  and  (17)  indicates  that  the  flow  velocity,  v,  is  in  a direction  of 
decreasing  values  of  z.  From  the  differentiation  of  Eq.(17)  it  can  be 

concluded  that  the  equation  is  an  exact  solution  for  a layer  of  infinite  depth. 
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However,  in  accordance  with  Terzaghi,  1942,  it  is  only  valid  for  the  time 
increment  indicated  if  the  depth  of  the  layer  is  limited.  Solutions  for  the 
transient  flow  condition  and  other  boundary  conditions  of  interest,  such  as 
cylindrical  and  spherical  surfaces  in  an  infinite  mass  can  be  found  in  Carslaw 
and  Jaeger,  1948,  pp  334-337,  and  a solution  for  a linearly  changing  surface 
pressure,  p = at,  where  a is  a constant,  was  presented  by  Clements,  1974. 

In  contrast  to  the  solutions  for  Eq  (11)  obtained  in  Eqs.(15)  through  (17), 
Equation  (12)  would  apply  if  the  pressure  difference  between  the  pore  pressures 
in  the  surrounding  soil  mass  and  the  interior  of  the  building  persists  for  a 
long  period  of  time,  so  that  pressure  gradients  can  stabilize.  Steady-state 
solutions  were  developed  for  many  boundary  conditions  using  the  semi-graphical 
technique  of  flow  nets  [i.e.  Cedergren,  1967,  Harr,  1962],  and  other  techniques, 
such  as  the  method  of  "fragments"  [Pavlovski,  1956,  Harr,  1962],  or  appropriate 
computer  programs.  Some  closed-form  solutions  were  also  developed  [i.e.  Dupuis, 
1863] . A solution  which  seems  appropriate  for  the  radon  entry  into  basements  was 
presented  by  DSMA  ADCON  Ltd.,  1983,  based  on  an  analogy  between  Ohm's  law.  The 
same  equation  is  also  used  by  Mowris  and  Fisk,  1988,  to  calculate  radon  entry 
rates  into  houses.  If  it  is  assumed  that  radon  entry  is  through  a linear  gap 
between  the  basement  slab  and  the  perimeter  wall,  then: 

^Lkp^ 

Q = — — ...(18) 

H cosh'^(D/a) 

where:  Q = rate  of  flow  (m^-s*^) 

Pg  = total  pressure  difference  (Pa) 
a = half  width  of  the  crack  (m) 

D = depth  below  surface  (m) 

L = length  of  crack  (m) 

DSMA  ADCON  make  the  case  that  for  D > 0.5  m and  a crack  width  of  10'^  m which  is 
typical  for  anticipated  conditions  this  equation  is  not  very  sensitive  to  depth 
change.  If  /i  is  taken  as  1.8*10’^  Pa-s,  and  cosh“^  (D/a)  as  > 8.3,  then: 

k L p^ 

Q < — __ — . . . . (19) 

4.8-10-5 
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Since  this  model  does  not  account  for  the  geometry  of  the  house  foundation  it  is 
suggested  by  DSMA  ADCON  to  add  1.5  m to  the  depth  of  the  crack  to  account  for 
the  foundation.  The  authors  justify  the  recommendation  by  a graphical 
construction  of  flow  lines. 

Steady-state  flow  is  also  assumed  in  a finite  difference  model  developed  by 
Oliveira-Loureiro , 1987.  This  model  considers  several  parameters  simultaneously 
by  solving  Eq(5)  in  a soil  block  extending  10  m from  the  outer  walls  of  the 
building.  A sensitivity  study  is  then  conducted  varying  one  parameter  at  a time 
while  the  other  parameters  are  kept  constant.  Permeability  is  one  of  the 
parameters.  As  expected,  it  was  concluded  that  if  soil  permeability  were  the 
only  factor  affecting  the  soil-gas  flow  there  would  be  a linear  correlation 
between  soil-gas  flow  and  permeability  k.  But  the  soil-gas  flow  and  the  radon 
flux  are  also  affected  by  other  variables  which  are  subsequently  discussed.  It 
must  be  noted,  however,  that  the  findings  of  the  study  have  no  experimental 
verification,  that  a steady-state  flow  model  was  used  which  represents  a lower 
bound  for  convective-flow  effects  and  that  many  other  simplifying  assumptions 
were  made.  Thus  the  conclusions  are  at  best  qualitative. 

First,  convective -flow  effects  are  limited  by  the  soil-gas  entry  mechanism.  The 
study  assumed  a crack  and  varied  the  crack  width  from  0.5  to  10  mm.  It  was 
found  that  only  for  cracks  wider  than  5 mm  was  there  no  significant  pressure 
drop  across  the  crack.  For  crack  widths  smaller  than  1 mm  the  crack  width  as 
well  as  k affected  the  flow  volume  for  k > 10’^^  m^ . Furthermore,  the  diffusive 
component  of  Eq(5)  dominated  the  radon  flux  for  k < 10"^^  m^ . For  soil 
permeabilities  k > 10”^^  to  10*^°  dilution  of  the  radon  concentration  in  the 
soil  gas  became  significant  and  was  very  significant  at  the  upper  limit  of  this 
permeability  range.  Thus  it  can  be  seen  that  even  with  a steady-state  flow 
model,  radon  flux  was  a linear  function  of  k only  for  a narrow  range  of 
permeabilities , 

The  following  overall  conclusion  on  the  effect  of  k on  the  convective  flux  of 
soil  gas  can  be  deduced  from  Eqs . 13,  through  19: 

(1)  For  transient  flow,  v = f(k°-^);  for  steady-state  flow  v = f(k). 
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(2)  At  t = ®,  and  z = 0,  v(pp,t)  is  a minimum.  Thus  the  steady-state  flow 

velocity  is  a lower  bound  at  the  point  where  the  initial  pressure  drop 

occurs . 

These  conclusions  must  be  considered  when  the  correlation  between  permeability 
and  the  rate  of  radon  entry  into  houses  is  assessed. 

3.  VARIABLES  AND  THEIR  MEASUREMENT 

3.1  Radon  Generation  in  Soil  Pores 

The  rate  -of  radon  generation  in  the  pore  space  can  be  measured  by  the  quantity 

^max  [Eq(4)].  The  variables  on  which  depends  are  Aj^^ , the  radium  activity 

concentration  in  the  soil  grains,  7^,  the  dry  density  of  the  soil,  and  f,  the 
emanation  coefficient.  These  variables  are  subsequently  discussed. 

Radium  Activity  Concentration 

The  quantity  Aj^^^  is  a measure  of  the  radium  content  of  the  soil.  Radium  (Ra^^®) 
has  a half-life  of  1600  years.  Thus,  for  the  time  frame  of  the  problem 
considered  the  radium  activity  concentration  is  an  invariable  quantity. 

A considerable  amount  of  data  on  radium  activity  concentrations  for  various 
soils  is  available  [i.e.  Nazaroff  et  a.,  1988,  DSMA  Atcon,  1983].  The  data 

suggest  that  the  radium  activity  concentration  in  surface  soils  ranges  typically 
from  10  to  100  Bq-kg“^.  However,  values  in  excess  of  1,700  Bq-Kg'^  have  been 
recorded  for  uranium  rich  soils.  A range  of  values  from  0.4  to  3,500  Bq-kg'^ 
were  recorded  for  rocks.  Values  tend  to  be  high  for  shales  and  schists 
[Gundersen  et  al.  , 1988]  and  tend  to  decrease  with  silica  content  in  igneous 

rocks  [Nazaroffet  al.,  1988]. 

The  Ra^^®  activity  concentration  can  be  measured  in  the  laboratory  by 
determining  the  Bi^^^  emission  of  samples  which  have  been  stored  long  enough  for 
the  Bi^^^  to  come  into  equilibrium  with  its  radium  precursor.  Gamma  spectroscopy 
is  used  to  make  the  measurement. 


9 


More  rapid  but  less  accurate  measurements  can  be  performed  in  the  field. 
Measurements  of  radium  activity  concentration  at  a shallow  depth  (less  than  0.5 
m)  can  be  made  by  detectors  suspended  above  the  ground  surface.  There  are  two 
problems  affecting  the  accuracy  of  these  measurements:  (1)  The  measured  emission 
decreases  with  increasing  moisture  contents.  This  error  could  be  corrected  by 
determining  moisture  contents.  (2)  Radium  concentrations  are  inferred  from 
Bi^^^  measurements.  However  radium  decays  directly  to  radon  gas,  which  in  turn 
decays  to  Po^^®,  a metal  which  rapidly  attaches  itself  to  a nearby  surface  where 
it  decays  to  Bi^^'* . Near  the  soil  surface,  part  of  the  radon  gas  escapes  to  the 
atmosphere.  Thus,  near  the  surface,  a measurement  of  Bi^^^  is  likely  to 
underestimate  the  radium  concentration  by  approximately  20%.  This  is  another 
error  which  needs  to  be  compensated  for,  probably  by  estimating  the  effect. 
Measurements  at  greater  depth  can  be  made  by  borehole  probes  which  have  been 
developed  for  uranium  detection. 

Gamma  ray  spectroscopy  is  carried  out  either  by  sodium  iodide  (Nal) 
scintillating  detectors  or  by  solid-state  diode  detectors  which  are  usually  made 
of  germanivim  (Ge)  . The  Nal  detector  measurements  have  lower  resolution  than  the 
Ge  detector  measurements,  but  Nal  detectors  operate  at  ambient  temperatures  and 
are  less  expensive  than  the  Ge  detectors.  Ge  detectors  operate  at  cryogenic 
temperatures.  While  the  precision  of  the  Ge  detector  measurements  is  much 
greater  than  that  of  the  Nal  measurements,  the  overall  measurement  accuracy  is 
not  necessarily  much  improved  because  of  the  previously  discussed  uncertainties 
associated  with  insitu  measurements. 

Dry  Density  of  Soils 

Dry  densities  of  soils  which  are  of  interest  in  this  study  range  from  1,200 
kg-m'®  for  loose  micaceous  sand  to  2,300  kg-m'®  for  silty  sands  with  gravels 
(mostly  glacial  tills).  Data  on  various  soil  types  are  given  in  Tables  la  and  b. 
The  data  tabulated  were  developed  from  different  information  bases.  The  95% 
modified  proctor  densities  are  for  ASTM  D1557  which  represents  compaction 
greater  than  that  normally  used  by  the  building  construction  industry  (ASTM 
D698)  and  should  be  considered  an  upper  bound  for  compacted  soils.  The  dry 
density  of  the  soil  in  situ  can  be  considered  an  invariable  quantity,  even 


10 


TABLE  la:  Typical  Soil  Properties 
(After  Lambe  and  Whitman,  1969) 


Soil  Type 

Void 

Ratio 

Porosity 

O 

Dry  Density  kg-m 

®max 

®min 

^max 

’^min 

Vdmin 

Vdmax 

Uniform  Spheres 

0.92 

0.35 

0.476 

0.26 

Standard  Ottawa  Sand 

0.80 

0.50 

0.44 

0.33 

1,470 

1,760 

Clean  Uniform  Sand 

1.00 

0.40 

0.50 

0.29 

1,330 

1,890 

Uniform  Inorganic  Silt 

1.10 

0.40 

0.52 

0.29 

1,280 

1,890 

Silty  Sand 

0.90 

0.30 

0.47 

0.23 

1,390 

2,030 

Fine  to  Coarse  Sand 

0.95 

0.20 

0.49 

0.17 

1,360 

2,210 

Micaceous  Sand 

1.20 

0.40 

0.55 

0.29 

1,220 

1,920 

Silty  Sand  and  Gravel 

0.85 

0.14 

0.46 

0.12 

1,430 

2,340 

TABLE  lb:  Typical  Dry  Densities  of  Soils 

(After  Holtz  and  Kovacs , 1981) 


Soil  Type 

0 

Dry  Density,  kg-m”-^ 

Vdmin 

Vdmax 

Sands  and  Gravels 

1,500 

2,300 

Silts  and  Clays 

600 

1,800 

Glacial  Tills 

1,700 

2,300 

Crushed  Rock 

1,500 

2,000 

Peats 

100 

300 

Organic  Silts  and  Clays 

500 

1,500 

Typical  Values  for  95%  Modified  Proctor  Density:  (after  Holtz  and  Kovacs,  1981) 


Clean  Uniform  Sand  (fine  to  medium) 

Well  Graded  Silty  Sand 

Well  Graded  Clean  Fine  to  Course  Sand 
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1,770  kg-m'3 

1,880 

2,030 


though  minor  changes  may  occur  as  a result  of  applied  surface  loads.  However  if 
the  soil  is  excavated  and  subsequently  replaced  as  fill  the  density  is  likely  to 
change.  Granular  soils  (gravels,  sands  and  silts)  compacted  in  accordance  with 
prevailing  standards  are  likely  to  have  densities  greater  than  their  original  in 
situ  density.  However  finegrained  soils  (silty  clays  and  clays)  and  some 
glacial  tills  may  have  densities  which  are  smaller  than  the  original  in  situ 
density. 

Several  methods  are  used  to  measure  dry  density.  For  cohesive  (fine  grained) 
soils  which  are  not  desiccated,  in  situ  density  as  well  as  water  content  can  be 
determined  in  the  laboratory  from  undisturbed  samples,  usually  taken  by  pressing 
a thin  wall  tube  into  the  soil.  However  undisturbed  samples  are  difficult  to 
take  in  granular  soils,  particularly  sands  and  gravels,  and  in  dry  cohesive 
soils,  and  density  determinations  from  disturbed  samples  (such  as  split  spoon 
samples  taken  in  the  Standard  Penetration  Test  (SPT)  [ASTM  D1586]  are  likely  to 
be  inaccurate.  Soil  density  can  be  measured  insitu  by  removing  a known  volume 
of  soil  and  determining  its  dry  weight  in  the  laboratory  [ASTM  D1556].  This 
method  is  somewhat  cumbersome  and  can  only  be  performed  on  the  soil  surface  or 
on  the  bottom  of  a test  pit  large  enough  for  a person  to  enter,  but  it  is 
commonly  used  in  compaction  control.  Nuclear  density  tests  [ASTM  D2922]  are  more 
rapid  and  less  cumbersome  than  the  ASTM  D1556  method.  The  nuclear  method  also 
usually  combines  density  and  moisture  content  measurements  (the  methods 
described  are  used  for  compaction  control  where  it  is  also  important  to 
determine  moisture  content).  Density  is  measured  by  7-radiation  scatter,  using 
radium  or  a radioactive  isotope  of  cesium  as  a source.  Moisture  content  is 
measured  by  neutron  scatter  caused  by  hydrogen  atoms.  Typically  americium- 
beryllium  isotopes  are  used  as  neutron  source.  Like  ASTM  D1556,  the  nuclear 
density  test  described  above  is  used  for  compaction  control  and  is  performed  on 
the  soil  surface.  However,  nuclear  borehole  probes  have  been  built  [i.e. 
Vaswani , 1975 ] . . 

It  should  be  noted  that  in  view  of  the  uncertainties  associated  with  estimates 
and  measurements  of  the  variables  associated  with  the  radon  generation  and 
transport  problem,  great  accuracy  in  the  measurement  of  the  dry  density  of  the 
soil  may  not  be  necessary  (however,  as  will  be  discussed  later,  accuracy  in  the 
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determination  of  the  moisture  content  is  important) . A reasonably  good  estimate 
can  be  made  with  the  aide  of  Table  1,  a definition  of  the  soil  type,  and  some 
information  from  penetration  tests,  such  as  the  SPT,  the  cone  penetration  test 
(CPT)  [ASTM  D3441)  or  the  helical  probe  test  (HPT)  [Yokel  and  Mayne , 1986).  The 
dry  density  of  man-made  fill  material  can  also  be  estimated  on  the  basis  of 
information  on  the  procedure  used  to  place  the  fill,  or  from  field  density  tests 
performed  during  the  placement  of  the  fill. 

Emanation  Coefficient  (f) 

While  the  radium  content  and  the  dry  density  of  the  soils  can  be  considered 
invariable  quantities  for  practical  purposes,  the  emanation  coefficient,  and 
thus  also  the  quantity  > strongly  depends  on  the  moisture  content  of  the 
soil.  The  mechanism  by  which  radon  atoms  are  emitted  to  the  pore  space  is 
discussed  by  Nazaroff  et  al.,  1988  and  Tanner,  1964  and  1980.  Since  linear 
momentum  is  conserved,  Rn^^^  atoms,  when  created  by  the  a decay  of  radium, 
posses  a kinetic  energy  of  86  keV.  The  atoms  travel  from  their  site  of 
generation  until  their  energy  is  transferred  to  the  host  material  for  an  average 
distance  which  depends  on  the  density  of  the  host  material.  This  distance  would 
be  0.02-0.07  /jm  in  common  minerals,  0.1  /xm  in  water,  and  63  jj-m  in  air.  Thus  on 
the  average  the  distance  at  which  the  path  of  the  radon  atom  terminates  in  the 
soil  gas  is  630  times  greater  than  the  average  travel  distance  in  water,  and  in 
many  instances  exceeds  the  travel  distance  available  in  the  pore  space.  Thus 
when  the  pore  space  is  filled  with  gas,  many  of  the  radon  atoms  emitted  traverse 
the  pore  space  and  embed  themselves  in  soil  grains.  Only  a relatively  small 
portion  of  the  radon  atoms  either  terminate  in  the  pore  space  or  return  to  the 
pore  space  by  diffusion.  If  the  pores  are  filled  with  water,  a much  larger 
portion  of  the  radon  atoms  emitted  remains  in  the  pore  space. 

Attempts  have  been  made  to  model  this  problem  mathematically,  considering  the 
three  components  of  the  emanation  coefficient,  direct  recoil,  indirect  recoil 
and  diffusion,  and  assuming  spherical  intact  soil  grains.  Interestingly,  it  was 
found  that  the  mathematical  model  yields  much  smaller  emanation  coefficients 
than  those  measured  in  the  laboratory  [Nazaroff  et  al.,  1988,  Tanner,  1980].  The 
reasons  advanced  for  the  discrepancy  are;  (1)  that  the  radium,  rather  than  being 
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distributed  uniformly  in  the  mass  of  the  soil  grains  is  concentrated  in 
secondary  crusts  near  the  grain  surface;  and  (2)  that  the  crystalline  structure 
in  the  vicinity  of  newly  formed  radon  atoms  is  damaged  by  chemical  corrosion  and 
radiation,  permitting  the  atoms  to  more  readily  migrate  into  the  pore  space  than 
the  density  of  the  soil  grains  would  otherwise  indicate.  Some  confirmation  for 
these  hypotheses  is  provided  by  measurement  of  emanation  coefficients  in  the 
Island  of  Hawaii  (Wilkening,  1974).  Emanation  coefficients  measured  on  samples 
from  lava  fields  of  recent  geological  origin  were  on  the  order  of  0.02,  while 
those  from  thin  organic  soils  were  on  the  order  of  0.55  and  those  from  deep 
agricultural  soils  0.7. 

Figure  1 (taken  from  Nazaroff  et  al.,  1988)  shows  a correlation  curve  between 
emanation  coefficients  and  moisture  content  developed  by  Strong  and  Levins,  1982 
for  a sample  of  uranium  mill  tailings.  The  curve  indicates  that  the  emanation 
coefficient  for  a dry  sample  was  approximately  22  % of  that  for  a saturated 
sample.  Initially  the  emanation  coefficient  increased  rapidly  with  the  moisture 
content  until  the  moisture  content  reached  approximately  25%  of  saturation,  at 
which  point  the  coefficient  is  80%  of  that  measured  at  saturation.  Some 
additional  data  are  tabulated  by  Nazaroff  et  al.,  1988,  p 74.  DSM  Atcon,  1983 
tabulated  data  derived  from  measurements  on  oven-dried  samples  (Table  4) . They 
conclude  that  typically  values  of  f for  oven-dried  soils  are  0.15  to  0.20  and 
recommend  to  assume  f values  for  moist  soils  to  be  "virtually  constant"  at  0.35. 
Unfortunately  there  is  only  a very  limited  data  base  to  justify  this  suggestion. 
For  instance  Greiner,  1985  measured  the  effect  of  moisture  content  on  the 
emanation  coefficient  from  ground  coal  and  concluded  that  the  emanation 
coefficient  for  a dry  sample  was  very  low. 

The  emanation  coefficient  can  be  measured  on  soil  samples  which  were  sealed  for 
long  enough  for  C to  reach  in  the  pore  space,  at  least  14  days  (3.7  half 
lives)  . A gas  volume  can  then  be  withdrawn  and  expanded  into  a Lucas  flask  to 
measure  the  radon  concentration.  An  efficient  measurement  method  for  for 
saturated  conditions  using  very  small  soil  samples  and  liquid  scintillation  was 
described  by  Wadach  and  Hess,  1985.  Both,  f for  dry  and  for  saturated 
conditions  are  invariant  quantities,  however  the  emanation  coefficient  in  the 
field  at  any  point  in  time  depends  on  the  moisture  content  of  the  soil.  The 
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Figure  1:  Effect  of  Moisture  Content  on  the  Emanation  Coefficient  for  a sample 
of  Uranium  Mill  Tailings  (from  Strong  and  Levins,  1982). 
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emanation  coefficient  for  the  dry  condition  is  a lower  limit.  The  maximum  radon 
flux  in  the  field  would  occur  in  a moist  condition,  when  there  is  enough 
moisture  content  to  substantially  increase  the  radon  generation  over  that  in  a 
dry  condition  but  not  enough  moisture  to  substantially  reduce  the  gas 
permeability  of  the  soil.  The  moisture  content  associated  with  this  condition 
depends  on  the  particle  size  distribution  in  the  soil  and  falls  somewhere 
between  the  field  capacity,  where  the  capillary  component  of  the  pore  water 
(pore  water  which  is  retained  by  surface  tension  and  does  not  drain  by  gravity 
forces)  is  retained  and  the  wilting  point,  where  only  the  hygroscopic  component 
of  the  pore  water  remains  in  the  pores  (electrostatically  adsorbed  pore  water 
which  is  retained  for  a long  period  of  time  in  otherwise  dry  soil) . The 
approximate  moisture  contents  associated  with  these  conditions  are  given  in 
Table  2 [after  Nazaroff  et  al.,  1988].  Except  in  arid  regions,  the  hygroscopic 
water  is  retained  most  of  the  time.  However  according  to  Nazaroff  et  al.,  1988, 
this  portion  of  the  pore  water  does  not  significantly  increase  the  emanation 
coefficient.  In  most  soils  the  capillary  water  is  retained  at  some  height  above 
the  groundwater  table. 

TABLE  2:  Typical  Water  Retention  Characteristics  of  Soils 
(taken  from  Nazaroff  et  al.,  1988) 


Soil  Type 

Grain  Size 

mm 

Porosity 

Saturation  (percent) 

Field  Cap. 

Wilting  Pt. 

Sand 

2.00  - 0.06 

0.4 

15 

5 

Silt 

0.06  - 0.002 

0.5 

58 

10 

Clay 

< 0.002 

0.6 

68 

48 

3,2  Radon  Transport  in  Soils 

The  variables  controlling  the  radon  flux  are  the  diffusion  coefficient  and  the 
gas  permeability  coefficient.  For  convenience,  the  formulations  in  this  report 
deal  with  bulk,  rather  than  effective  coefficients.  Thus  while  the  actual  radon 
flux  takes  place  within  the  part  of  any  cross  section  area  which  is  occupied  by 
the  pores,  the  bulk  coefficients  are  averaged  over  the  entire  cross  section 
area.  The  diffusion  and  permeability  coefficients  are  subsequently  discussed. 
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Diffusion  Coefficient 


While  radon  exhalation  from  the  exposed  soil  surface  to  the  atmosphere  is  to  a 
large  extent  the  result  of  diffusion,  the  dominant  transport  mechanism 
associated  with  the  entry  of  radon  into  buildings  is  convection,  unless  the  gas 
permeability  coefficient  is  very  low,  a situation  which  is  usually  associated 
with  low  radon  entry  rates.  Since  the  diffusion  coefficient  varies  primarily 
with  moisture  content,  but  does  not  vary  much  with  soil  type,  it  is  not 
recommended  in  this  report  to  measure  the  diffusion  coefficient.  For  practical 
purposes  it  is  considered  adequate  to  use  Eq(9)  for  estimating  the  effective 
diffusion  coefficient.  To  obtain  the  bulk  diffusion  coefficient  the  results  of 
Eq(9)  must  be  multiplied  by  n.  An  upper  bound  for  the  effective  diffusion 
coefficient  is  the  diffusion  coefficient  of  radon  in  open  air,  , which  is 
1.2- 10”^  m^-s'^.  A lower  bound  in  saturated  soil  is  2*  10'^°  m^-s'^.  Note  that  as 
a result  of  diffusion  the  radon  concentration  near  the  soil  surface  is  less  than 
that  at  greater  depths.  In  an  ideal  radon  concentration  profile  where  C = 0 at 
the  soil  surface  and  C = at  an  infinite  depth  C would  be  0.63  at  the 
depth  of  li  which  is  approximately  1 m. 

Permeability  Coefficient 

The  permeability  coefficient  is  a measure  of  the  frictional  resistance  to  the 
flow  of  soil  gas  in  the  pores.  It  depends  on  the  shape  and  size  of  the  pores. 
Since  in  most  soils  the  pores  are  at  least  partially  filled  with  water,  the  gas 
permeability  is  strongly  dependent  on  the  moisture  content  of  the  soil.  The  gas 
permeability  of  dry  soils  and  the  water  permeability  of  saturated  soils  can  be 
considered  invariant. 

Since  in  the  mathematical  formulations  previously  presented  the  soil  was  assumed 
isotropic  with  respect  to  permeability,  k could  be  treated  as  a scalar  quantity. 
However,  because  of  the  manner  in  which  soils  in  nature  (and  also  during 
construction)  are  deposited  and  pre-consolidated,  the  horizontal  permeability  is 
in  many  instances  greater  than  the  vertical  permeability.  The  ratio  between  the 
horizontal  and  vertical  permeabilities  ranges  from  about  three  for  soils  which 
are  pre-consolidated  but  not  stratified  to  more  than  an  order  of  magnitude  and 
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sometimes  two  orders  of  magnitude  for  stratified  sedimentary  soils  such  as 
varved  clays.  In  some  residual  soils,  which  formed  by  the  decomposition  of  the 
local  rock  formation  rather  than  by  deposition  of  particles  transported  from 
other  locations  and  which  were  not  heavily  pre-consolidated  during  their 
geological  history,  the  vertical  permeability  exceeds  the  horizontal 
permeability.  Note  that  the  northern  part  of  the  United  States  was  glaciated 
during  its  recent  geological  history  and  that  therefore  most  soils  in  these 
regions  were  pre  consolidated  by  thick  sheets  of  ice, 

A considerable  amount  of  data  is  available  for  saturated  water  permeability  from 
field  and  laboratory  tests.  These  data  can  be  used  to  approximately  estimate 
dry  gas  permeability.  However,  it  is  important  to  recognize  that  for  flow  of 
water  through  soils  the  boundary  velocity  at  the  surface  of  the  soil  grains  is 
zero,  whereas  for  gases  it  is  not  zero  [Corey,  1988,  p.98].  This  phenomenon  was 
first  observed  by  Klinkenberg,  1941,  and  is  referred  to  as  the  "Klinkenberg 
effect"  or  slipflow.  The  value  of  k derived  from  water  permeability  is  a lower 
bound  which  is  approached  by  gas  flows  at  very  large  pressures.  At  atmospheric 
pressures,  the  gas  permeability  will  be  approximately  equal  to  the  water 
permeability  for  coarse  sands.  It  could  be  as  large  as  twice  the  water 
permeability  for  particle  sizes  in  the  silt  range  [Corey,  1988],  In  view  of  the 
fact  that  permeabilities  range  over  several  orders  of  magnitude  and  that  for 
finegrained  soils  the  radon  flux  is  primarily  attributable  to  diffusion  or 
convection  through  fissures,  the  increase  in  permeability  attributable  to  the 
Klinkenberg  effect  is  not  very  significant.  Nevertheless,  gas  permeabilities 
derived  from  data  on  water  permeability  should  be  considered  lower  bound  values. 


Since  Darcy's  coefficient  for  water  permeability  is  given  for  the  gradient  of 
the  water  head  and  has  the  dimension  of  L-T”^  (velocity),  it  can  be  converted 
to  k as  follows: 


, dH  kyw  dH 

V = kw  Xo  - 

^ dS. 


di 


k > kw  — 
7w 


where  yw,  /iw  and  kw  are  the  unit  weight,  dynamic  viscosity  and  permeability 
coefficient  of  water. 


(a  18°C:  /jw  = 1.053-10'3  Pa-s  and  yw  = 9.8-103  N-m’^ 
therefore:  k (m^)  > kw  (m- s'^)  • 1 .07  • 10“^  (ra-s)  ...(20) 
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Some  data  derived  from  information  on  water  permeability  are  given  in  Table  3 
which  describes  soil  characteristics  and  some  methods  of  determining  water 
permeability.  The  data  were  originally  presented  by  Casagrande  and  Fadum,  1944. 
The  laboratory  methods  listed  are  constant  head  permeameter,  falling  head 
permeameter,  horizontal  capillarity  test  and  consolidation  test  (Lambe,  1958). 
Most  of  the  laboratory  tests  can  be  used  to  determine  the  water  permeability  of 
saturated  samples  (though  in  some  instances  it  may  be  difficult  to  drive  all  the 
soil  gas  from  the  pore  space) . However  the  horizontal  capillarity  test  measures 
the  permeability  of  the  capillary  portion  of  the  pore  space  (at  approximately 
80%  saturation)  which  is  approximately  33  to  50%  of  the  saturated  permeability 
[Lambe,  1958].  The  permeability  of  the  pores  larger  than  those  in  the  capillary 
range  would  be  the  difference  between  the  permeability  at  saturation  and  the 
permeability  obtained  in  the  horizontal  capillarity  test.  This  value  is  of 
interest,  because  the  emanation  coefficient  of  partially  saturated  soils  is  much 
higher  than  that  of  dry  soils  and  the  degree  of  saturation  in  the  horizontal 
permeability  test  may  be  similar  to  that  below  basement  foundations  which  are 
close  to  the  phreatic  surface,  where  the  soil  is  partially  saturated  as  a result 
of  capillary  suction. 

The  only  field  permeability  test  listed  in  Table  3 is  the  pumping  test,  which 
can  only  be  conducted  below  the  groundwater  table.  While  this  latter  test  is 
not  practical  for  the  purpose  site  classification  for  radon,  it  provided  a 
large  volume  of  reliable  empirical  permeability  data.  Many  other  types  of  in 
situ  water  permeability  tests,  mostly  borehole  tests  have  been  developed.  An 
overview  and  assessment  of  available  test  methods  is  provided  by  Milligan,  1975. 
The  permeability  of  soils  can  also  be  measured  electrically  and  characterized  by 
a "formation  factor"  which  is  the  ratio  of  the  conductivity  of  the  electrolyte 
occupying  the  pore  space  of  a saturated  soil  sample  to  the  conductivity  of  the 
soil-solution  mix.  This  measurement  can  be  performed  in  situ  and  can  also  be 
used  to  determine  anisotropy.  The  problem  with  this  measurement  is  that  the 
soil  must  be  saturated  [Arulmoli  et  al.,  1985]. 
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Table  3.  Soil  Permeability  and  Drainage  Characteristics 
(based  on  data  presented  by  Holtz  and  Kovacs,  1981,  Peck  et  al.,  1974) 
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A test  of  considerable  practical  interest  is  a percolation  test  developed  by  the 
Federal  Housing  administration  (FHA)  for  the  purpose  of  designing  subsurface 
wastewater  absorption  fields  [HUD  Minimum  Property  Standards,  1973],  The  test 
is  required  by  most  jurisdictions  in  areas  not  served  by  municipal  sewage 
disposal  systems.  It  is  performed  in  an  open  auger  hole  or  pit.  The 
percolation  rate  is  measured  over  a four-hour  period  after  the  soil  was 
saturated  for  a specified  period  of  time.  While  there  are  presently  no  data 
correlating  percolation  rate  with  water  permeability,  such  data  can  be 
generated.  Considerable  economy  could  be  realized  by  utilizing  these  test 
results  for  radon  risk  evaluation. 

The  gas  permeability  values  listed  in  Table  3 were  converted  from  water 
permeability  using  Eq  (20)  and  therefore  represent  lower-bound  values.  The 
vertical  double  lines  represent  divisions  between  easily  recognizable  soil 
categories  (gravels;  sands;  silts  and  clays)  and  are  considered  "benchmark" 
values . 

Of  particular  interest  are  the  modified  "impervious"  soils.  The  permeability  of 
these  soils  cannot  be  correlated  with  their  grain  size  distribution,  because 
their  permeability  is  caused  by  structural  defects  and  exceeds  their  intact 
permeability  by  up  to  several  orders  of  magnitude.  The  range  of  values  given  in 
Table  3 extends  over  four  orders  of  magnitude  and  is  derived  from  empirical 
data.  The  most  common  causes  of  structural  defects  are  cycles  of  swelling  and 
shrinkage,  shrinkage  caused  by  desiccation,  and  effects  of  vegetation.  In  the 
first  category  are  the  swelling  clays  in  semi-arid  regions  in  the  U.S.  These 
soil  deposits  shrink  and  swell  seasonally  and  develop  fissures  in  the  dry 
season.  As  a result  of  the  cyclic  volume  changes  these  soils  are 
overconsolidated  and  develop  structural  defects  when  lateral  confinement  is 
removed  by  excavation.  Desiccated  clays  may  be  fissured  before  construction 
starts.  In  other  instances  clay  deposits  may  have  a high  enough  natural 
moisture  content  to  prevent  shrinkage  cracking,  but  will  eventually  desiccate 
when  covered  by  building  foundations.  Vegetation  effects  are  of  particular 
interest  in  slab-on-grade  buildings  which  rest  on  very  shallow  soil  deposits. 
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Shrinkage  cracking  occurs  locally  when  the  capillary  pressures  exceed  the 
tensile  strength  (cohesion)  of  the  soil.  The  propensity  of  clays  for  developing 
these  structural  defects  is  related  to  measurable  index  properties.  Soils 
identified  in  the  Table  4 as  having  a high  or  very  high  degree  of  expansion 
should  be  considered  susceptible  to  shrinkage  cracking.  Note  that  the  Atterberg 
limits  alone  can  provide  the  necessary  information. 

Table  4:  Estimate  of  Expansion  from  Classification  Tests 
(After  U.S.  Bureau  of  Reclamation,  1974) 


Degree  of 
Expansion 

Probable  Exp.(% 
of  volume) , dry 
to  saturated  ■* 

Colloidal 
Content 
(%  - lyum) 

Plasticity 
Index  (PI)** 

Shrinkage 
Limit  (SL)** 

Very  High 

> 30% 

> 28% 

> 35 

< 11 

High 

20  - 30% 

20  - 31% 

25  - 41 

7 - 12 

Medium 

10  - 20% 

13  - 23% 

15  - 28 

10  - 16 

Low 

< 10% 

< 15% 

< 18 

> 15 

* Expansion  is  under  a surcharge  of  6.89  Pa  (1  psi) 

**  Atterberg  Limits  according  to  ASTM  D4318 

A simple  test  that  could  be  used  to  identify  fissured  clays  is  the  "drying  test" 
[Yokel  and  Tucker,  1980].  A small  clay  sample  is  allowed  to  dry  and  its  dry 
strength  and  integrity  is  examined  manually. 

The  permeability  caused  by  structural  defects  in  soils  can  be  measured  in 
particular  instances,  but  cannot  be  quantified  in  terms  of  other  measurable  soil 
index  properties.  However  clays  susceptible  to  shrinkage  cracking  can  be 
identified  by  their  swelling  potential.  Studies  of  this  problem  have  been 
conducted  by  Herzog  and  Morse,  1984,  Long  et  al . 1982,  McKinley  et  al.  1978, 
Snow,  1969,  and  Bosscher  et  al.  1988. 

Gas  permeability  can  be  measured  by  in  situ  tests.  Methods  have  been  developed 
to  combine  these  tests  with  measurements  of  radon  concentrations  in  the  soil  gas 
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Figure  2:  Relative  Air  Permeability  of  a Loamy  Sand  (from  Nazaroff  and  Nero, 
1988) 
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[Tanner,  1987,  Kunz , 1988,  Sextro,  1988,  Surback  and  Filler,  1988].  While  these 
measurements  provide  valuable  information,  they  are  sensitive  to  the  moisture 
content  at  the  time  of  measurement.  Figure  2,  which  was  taken  from  Nazaroff  et 
al.,  1988  shows  correlations  between  moisture  content  and  gas  permeability  for  a 
loamy  sand  (silty  sand)  . It  can  be  seen  that  for  this  soil  a large  correction 
would  be  required  for  saturations  in  excess  of  55%,  and  that  for  smaller 
saturations  the  permeability  is  very  sensitive  to  the  moisture  content. 

Typical  values  for  gas  permeability  presented  by  Nazaroff  et  al.,  1988  are 
shown  in  Table  5.  The  information  in  the  table  is  in  part  derived  from  the  same 
sources  as  that  in  Table  3 and  thus  is  compatible  with  Table  3.  Some 
mathematical  models  were  developed  which  correlate  particle  size  distribution 
with  permeability  (refer  to  Nazaroff  et  al.,  1988,  61-65,  Leonards,  1962,  119- 
124)  . Two  models  are  discussed  because  they  can  be  used  for  site 
characterization.  Hazen,  1892,  derived  an  empirical  formula  on  the  basis  of 
data  for  sharp,  uniform  sands  used  for  filters.  His  classical  formula  considered 
particle  size  and  had  a correction  for  the  temperature  effects  on  viscosity. 
The  following  expression  can  be  derived  from  Hazen' s formula: 

k = cdio^  ...(21) 

where:  c = a coefficient  accounting  for  density 

dj^Q  = the  particle  size  which  is  equal  or  larger  than  10%  of  the 

particles  by  weight,  also  called  the  "effective"  particle  size. 

Hazen  found  that  his  formula  is  valid  for  sands  whose  uniformity  coefficient: 

*^60 

Cu  < 5 , where  Cu  = ^ — = the  uniformity  coefficient  which  is  a measure  of 

^10 

grain  size  variation.  Hazen  observed  that  the  highest  value  of  c was 
approximately  twice  the  lowest  value.  For  gas  permeability,  the  Hazen  equation 
can  be  written  as: 

k (m^)  = (0.7  to  1 . 4)  • lO'^dj^g^  where  d is  in  m ...(22) 

A good  average  value  for  the  bracketed  term  in  Eq(22)  is  1 . In  present 
engineering  practice  a range  of  values  from  1 to  1.5  is  often  suggested  [Peck  et 
al.,  1973],  rather  than  the  range  in  Eq(22)  which  was  originally  determined  by 
Hazen.  It  was  suggested  that  the  Hazen  formula  gives  only  an  indication  of  the 
order  of  magnitude  of  permeability  [Sowers  and  Sowers,  1970]  and  its  validity 
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Table  5:  Typical  Gas  Permeability  Values 

(From  Nazaroff  et  al.,  1988) 


Log  k (m^) 


- -7 


Clean  gravel 


-8 


c 

1 

e 

a 

n 

s a 

n 

d 

s , 

c 

1 

e 

a 

n 

s a 

n 

d 

and 

s 

r 

a 

V 

e 

1 m 

i 

X 

t u r e 

- -11 


^ Uniform  coarse  sand 


Uniform  medium  sand 


Clean  well  graded  sand  & gravel 
•<-  Uniform  fine  sand 


Very  fine  sands, 

Organic  and  inorganic 
silts , 

Mixtures  of  sand,  silts 
and  clay,  etc. 


-12 


-13 


-14 


Well  graded  silty  sand  and  gravel 


■*-  Silty  sand 
Uniform  silt 


Sandy  clay 


- -15 


Silty  clay 


Homogeneous  clays 


-16 


^ Clay 
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was  thought  to  be  limited  to  loose  clean  uniform  sand  (with  less  than  5%  of 
particles  passing  the  # 200  sieve)  [Leonards,  1962,  Holtz  and  Kovacs,  1981]. 
However,  data  presented  by  Milligan,  1976,  p,13  and  Fig. 6 indicate  very  good 
agreement  between  the  Hazen  formula  and  high  quality  laboratory  and  pumping  test 
data,  even  for  uniformity  coefficients  as  high  as  20.  As  Table  3 indicates,  the 
Hazen  equation  applies  in  the  range  of  k > 10‘^^m^  which  is  the  range  of 
permeabilities  where  convective  flow  tends  to  dominate  radon  transport. 


The  Hazen  formula  does  not  have  a density  term  (even  though  density  effects  are 
implicit  in  the  variation  of  the  coefficient  c)  . Taylor,  1948,  suggested  the 
following  correlation  between  the  void  ratio  e = n/(l-n)  and  the  permeability: 

ki  C^ei^  0362^ 

ko  1 + e,  ’ 1 + e. 


..(23) 


For  sands  Cj^  = C2 . Thus  the  two  coefficients  can  be  factored  out  of  the 
equation.  A simpler  relationship  is  frequently  used  [Sowers  and  Sowers,  1970]: 


k oc  e2  = [n/(l-n)  ]2  . . . (24) 

A more  sophisticated  mathematical  model  for  k was  presented  by  Rogers  and 
Nielson,  1988: 

k=n-(l+n)-|^  I j%xp(-1.4m-2m2)  . { 0 . 0034  • exp  [ 0 . 264  • S^  • ^3 -0 . 6m- ( 1 -0  . IS)  ] ) . . . (25) 

where:  m = moisture  saturation  (as  a fraction  of  full  saturation) 

d = geometric  mean  particle  diameter  (by  size) 

S = geometric  standard  deviation  of  particle  diameter  distribution 


The  equation  assumes  a log-normal  particle  size  distribution  (by  diameter) . This 
is  a realistic  assumption  for  many  natural  granular  materials.  The  first  half 
of  Eq(25)  is  for  a geometric  standard  deviation  of  10,  and  the  second  half  of 
the  equation  (inside  the  braces)  adjusts  the  calculated  values  for  other 
geometric  standard  deviations.  An  interesting  aspect  of  Eq(25)  is  that  d^Q 
represents  a close  approximation  to  the  geometric  mean  particle  diameter,  and 
is  approximately  the  geometric  standard  deviation.  On  this  basis  a comparison 
is  made  between  Eqs . (22)  and  (25),  using  S=3  and  5,  n = 0.5  and  m = 0.  The 
equation  is  dimensionally  consistent.  Thus,  if  d is  in  m,  k is  in  m^ . 
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It  can  be  seen  that  the  calculated  permeabilities  are  in  the  same  range,  even 
though  the  value  calculated  by  Eq(25)  for  S = 5 is  higher  than  the  upper  limit 
suggested  for  Eq  (20).  This  is  encouraging,  since  Eq(25)  was  derived  from  an 
entirely  different  data  base  than  that  used  for  Eq(22) . Eq(25)  requires 
consideration  of  porosity  and  moisture  content. 


The  term  for  moisture  content  in  Eq(25)  is  of  interest  because  it  could 
potentially  be  used  to  correct  field  measurements  of  gas  permeability  for  the 
moisture  content  at  the  time  of  the  measurement.  Eq(25)  can  be  rewritten  as: 


k=0.0034n(l+n)  |j “ exp (0 . 264 • S^- 33 ) { exp [ - 1 . 4m-2m5-0 . 6m( 1 -0 . IS) ] } ... (26) 


The  term  inside  the  braces  in  Eq(26)  is  a correction  for  the  moisture  effect. 
Discussions  with  the  authors  of  Ref [45]  indicate  that  the  moisture  term  is  based 
on  a limited  amount  of  laboratory  data  for  soils  ranging  from  sands  to  clays, 
together  with  an  analytical  model  for  gas  permeability.  Moisture  content 
effects  calculated  from  Eq(26)  are  plotted  in  figure  3,  together  with  the  data 
from  figure  2 and  a typical  curve  shown  by  Corey,  1988.  The  data  from  figure  2 
and  Corey's  curve  are  not  inconsistent,  because  Corey's  curve,  which  is  a plot 
of  Eq  3.34,i  in  Corey,  1988,  depends  on  the  "residual  saturation"  and  on  an 
exponential  coefficient  which  are  soil-specific  (the  data  in  figure  2 are 
associated  with  a higher  residual  moisture  content  than  Corey's  curve  shown  in 
Figure  3).  Figure  3 indicates  that  Eq  (24)  does  not  provide  a realistic 
correction  factor  for  the  effect  of  saturation.  More  data  are  needed  to  develop 
correction  factors  for  moisture  content  which  can  be  correlated  with  readily 
measurable  soil  index  properties. 


27 


Figure  3:  Comparison  of  data  and  predictions  of  the  effect  of  the  degree  of 
saturation  on  the  gas  permeability  of  soils. 
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3 . 3 Summary 


The  radium  activity  concentration  can  be  measured  in  the  laboratory  or  in  the 
field  using  7-spectroscopy.  It  varies  typically  from  10  to  100  Bq-kg"^  of  dry 
soil , 

The  dry  density  of  the  soil  can  be  measured  in  the  field  by  removing  a known 
volume  of  soil  and  measuring  its  dry  weight  in  the  laboratory  or  by  using  a 
nuclear  density  gage.  It  varies  typically  from  1,200  to  2,300  kg-m'^  and  can  be 
estimated  from  data  on  soil  classification  and  penetration  tests. 

The  emanation  coefficient  is  sensitive  to  the  moisture  content  of  the  soil.  It 
can  be  measured  in  the  laboratory,  but  it  has  been  suggested  [DSMA  Atcon,  1983] 
that  estimates  should  be  used.  It  typically  is  on  the  order  of  15  to  20%  for  dry 
soils  and  35%  for  wet  soils.  However,  in  some  instances  emanation  coefficients 
higher  than  70%  have  been  recorded. 

The  diffusion  coefficient  decreases  with  moisture  content.  It  can  be  estimated 
closely  by  Eq  (9)  which  was  developed  from  a large  data  base.  It  varies  from 
7-10~®  m^-s'^  for  dry  soils  to  2- 10"^°  m^-s”^  for  saturated  soils. 

The  permeability  coefficient  decreases  with  increasing  moisture  content.  The 
coefficient  of  gas  permeability  for  dry  soils  is  an  invariable  quantity,  it  can 
be  estimated  from  the  particle  size  distribution  and  dry  density  of  the  soil  or 
measured  in  the  field  or  in  the  laboratory  by  various  tests.  It  varies  from  10'^ 
m^  for  clean  gravels  to  10"^®  m^  for  clays. 

4.  RADON  SOURCE  POTENTIAL  INDEX 

4.1  Radon  Entry  By  Convective  Flow 

Radon  could  enter  a building  in  various  ways,  depending  on  the  type  of 
construction  and  the  operation  of  the  HVAC  system.  In  this  section  an 
assessment  of  the  relative  effect  of  the  subsoil  permeability  is  made. 
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Radon  entry  is  in  most  instances  through  openings  in  basement  slabs  and  walls. 
Typically,  basement  slabs  are  underlain  by  a sub-base  of  4 inches  or  more  of 
granular  material  which  has  a high  gas  permeability.  If  entry  is  through  cracks 
and  other  openings  in  the  basement  slab  or  through  joints  between  the  slab  and 
the  foundation  wall,  which  in  most  instances  would  be  the  case,  the  granular 
sub-base  serves  as  a reservoir  in  which  equipotential  conditions  are  assumed  in 
order  to  simplify  the  problem  (actually  some  pressure  gradients  have  been 
observed  [DSMA  Atcon,  1983]).  There  may  or  may  not  be  a significant  pressure 
drop  between  the  sub-base  and  the  interior  of  the  basement,  depending  on  the 
width  of  the  openings  through  which  the  soil  gas  enters  and  the  quantity  of  gas 
that  enters  through  the  openings . Radon  entry  could  also  occur  through  cracks 
in  basement  walls,  which  would  result  in  a different  convective -flow  pattern. 

The  quantity  of  soil  gas  that  enters  a house  as  a result  of  a pressure  gradient 
between  the  sub-base  and  its  supporting  and  surrounding  soil  depends  on  the 
convective  flow  generated  by  the  pressure  gradient,  which  in  turn  depends  on  the 
permeability  and  thickness  of  the  supporting  soil. 

As  noted  in  the  conclusions  drawn  from  theoretical  derivations  in  Section  2.2, 
the  quantity  of  soil  gas  supplied  in  a given  length  of  time  as  a result  of  a 
pressure  gradient  from  the  subsoil  to  the  house  is  a function  of  k°-^  under 
transient  flow  conditions  and  a function  of  k under  steady-state  flow 
conditions.  This  consideration  directly  affects  the  dimension  of  a radon  index 
containing  k.  It  was  also  noted  that  the  steady-state  flow  condition,  which 
must  follow  a period  of  transient  conditions,  is  a lower  limit.  Thus  it  is 
important  to  establish  an  approximate  time  scale  which  can  be  compared  with  the 
length  of  time  of  predicted  cycles  of  pressure  drop. 

If  it  is  assumed  that:  the  soil  permeability  is  k;  the  openings  in  the  basement 
slab  are  large  enough  not  to  generate  a significant  pressure  drop  between  the 
interior  of  the  building  and  the  sub-base  of  the  slab;  the  thickness  of  the 
supporting  soil  layer  is  H and  the  soil  rests  on  impermeable  bedrock;  the 
lateral  dimension  of  the  house  is  much  larger  than  the  thickness  of  the  soil 
layer;  and  the  pressure  drop  occurs  at  time  t^  and  continues  for  a length  of 
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time  t;  then  the  following  equation  can  be  derived  from  Eq(17)  for  1 -dimensional 
flow: 

t.  = 

where:  t^  = time  required  to  displace  fraction  a of  the  total  quantity  of  gas 

that  would  be  displaced  from  a prismatic  soil  column  of  height  H 
by  a pressure  drop  p^  below  the  initial  atmospheric  pressure  . 

Equation  (27)  is  independent  of  p^ . This  is  only  true  when  Darcy's  law  is  valid 
and  the  pressure  drop  is  much  smaller  than  the  atmospheric  pressure.  Figure  4 
shows  a plot  of  tg  values  as  a function  of  permeability,  calculated  for  an  a 
value  of  0.5,  a value  of  1.013-10^  Pa,  a fi  value  of  1.827- 10'^  Pa-s  and  an  n 
value  of  0.5. 

Note  that  it  would  take  about  five  hours  to  displace  half  the  pore  gas  that 
would  be  displaced  from  a 10-m  (30-ft)  thick  column  of  silty  sand  by  a given 
pressure  drop,  but  only  about  three  minutes  to  displace  the  same  fraction  of  the 
soil  gas  from  a 10-m  thick  column  of  well  graded  clean  sand  and  gravel.  However, 
it  would  take  only  one  tenth  of  this  time  to  displace  half  the  soil  gas  from  a 
3-m  (10-ft)  thick  layer,  which  would  be  about  half  an  hour  for  a 3-m  column  of 
silty  sand  and  about  20  seconds  for  a 3-m  colximn  of  well  graded  sand  and  gravel. 
The  time  scale  in  figure  4 only  applies  when  the  lateral  dimensions  of  the  house 
are  much  larger  than  the  height  of  the  soil  column.  In  reality,  the  effect  of 
horizontal  inflow  of  gas  into  the  soil  along  the  perimeter  of  the  house  would  be 
to  substantially  prolong  the  duration  of  the  transient  flow  conditions  for  the 
thicker  soil  layers. 

Some  data  on  observed  fluctuations  of  pressure  differences  between  the  inside 
and  outside  of  buildings  are  given  by  Turk  et  al.,  1988,  and  Hubbard  et  al., 
1988.  From  their  data  plots  it  appears  that  in  summer,  when  pressure 
differences  caused  by  diurnal  temperature  cycles  and  wind  conditions  fluctuate 
between  positive  and  negative,  transient  flow  conditions  would  predominate. 
This  is  also  indicated  by  the  results  of  mathematical  simulations  of  wind 
effects  by  Scott  Assoc.,  1985  in  which  it  was  found  (by  a transient  convective- 
flow  model)  that  radon  entry  is  proportional  to  k°-^^.  In  winter,  however, 
because  of  the  stack  effect,  there  was  an  average  negative  pressure  gradient 
from  the  outside  to  the  interior  of  the  building,  in  addition  to  pressure 
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PERMEABILITY  (m  2) 


TIME  (min) 


Figure  4;  Time  required  to  displace  by  1 -dimensional  flow  50%  of  the  soil  gas 
that  would  be  displaced  by  a pressure  drop  at  time  t = 0. 
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fluctuations  which  were  of  the  same  order  of  magnitude  as  the  average.  While 
there  is  a long-term  negative  pressure  gradient  into  the  house,  the  time  scale 
of  the  soil  response  indicates  that  there  will  always  be  pressure  fluctuations 
which  will  result  in  a significant  transient  flow  component. 

The  steady-state  flow  occurs  after  a time  t which  is  long  enough  for  the  flow  to 
stabilize.  While  the  transient  flow  in  its  initial  stages  would  be  primarily 
along  a flow  path  from  the  soil  underlying  the  sub-base  vertically  upward  into 
the  sub-base,  the  steady-state  flow  would  be  along  a path  from  the  surface  of 
the  soil  surrounding  the  building  toward  the  sub-base  through  its  surrounding 
and  underlying  soil.  This  latter  flowpath  would  substantially  contribute  to  the 
dilution  of  the  radon  concentration  in  the  soil  gas  transported  and  thereby  tend 
to  reduce  the  effect  of  k on  the  radon  source  potential. 

4.2  Expressions  for  Radon  Source  Potential. 

Several  diagnostic  procedures  for  the  evaluation  of  radon  source  potential  on 
the  basis  of  measured  soil  properties  have  been  proposed.  These  range  from 
simple  expressions  which  only  consider  invariant  index  properties  of  the  soil, 
to  more  complex  expressions  which  also  consider  the  characteristics  of  the 
building.  Most  procedures  utilize  in  situ  measurements  of  soil  permeability  and 
soil  gas  radon  concentration  which  are  sensitive  to  barometric  pressure, 
temperature,  and  moisture  conditions  at  the  time,  as  well  as  prior  to  the  time 
the  measurements  are  made . 


Eaton  and  Scott,  1984,  proposed  a radon  index  number: 


hAi 


•Ra 


RIN  = 


...(28) 


where : 


-logk 

h = air  change  period,  hours. 


In  a later  study  with  a mathematical  model  [DSMA  Atcon,  1985]  the  same  authors 
found  that  radon  concentration  varied  as  k°-^^  and  suggested  the  following 
index : 

R = iog(AR,)  + 0.45-log(k)  ...(29) 
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The  authors  suggest  that  if  is  expressed  in  pCi-g”^  and  k in  cm^  a number  of 
R = -2.6  would  result  in  average  radon  concentrations  of  2-3  pCi/i.  Larger 

numbers  would  indicate  higher  average  concentrations. 

Another  index  number  is  proposed  by  Kunz  et  al.,  1988  and  backed  up  with  data 

from  glacial  terrace  deposits  in  New  York  State: 

RIN  = 10-Cm-km°-^  . . . (30) 

where:  Cm  = measured  radon  concentration  in  soil  gas  (pCi-i”^) 

km  = measured  gas  permeability  in  cm^ 

The  equation  is  calibrated  to  give  average  indoor  radon  concentrations  in  pCi/i 
and  has  been  calibrated  to  match  field  data.  For  the  data  examined  the  equation 
seems  to  provide  a reasonably  good  prediction  of  average  measured  indoor  radon 
concentrations.  Kunz  et  al.  also  report  several  interesting  observations.  An 
inverse  relationship  between  the  measured  radon  concentrations  in  the  soil  gas 
and  the  measured  permeabilities  was  observed.  This  could  be  attributable  to 
differences  in  moisture  content  (a  higher  moisture  retention  in  soils  with  more 
fines)  or  to  the  effect  of  pressure-gradient  induced  convective  flow  from  the 
atmosphere  into  the  soil.  It  was  also  established  that  the  radon  concentrations 
in  the  soil  gas  increased  with  the  distance  from  the  house,  probably  because  of 
dilution  resulting  from  convective  flow.  Both  of  these  observations  indicate 
that  the  effects  of  increased  convective  flow  are  partially  offset  by  decreased 
radon  concentration  in  the  soil  gas. 

For  the  same  soil  data,  an  index  number  can  be  derived  in  terms  of  invariant 
soil  index  properties.  Assuming  that: 

k = 2km  (a  moisture  correction;  "m"  means  "measured") 

^max  “ ~ 0.35  ( b came  out  to  be  approximately  2.4) 

And  that  for  the  conditions  examined  = 1,800  kg-m'^  and  n = 0.4,  the 
following  approximate  radon  index  can  be  derived  from  the  data  presented  by 
Kunz  et  al . : 

Y = 1.157dAR3-n-i-k°-5  ...(31) 

where:  is  in  kg-m"^;  is  in  Bq-kg"^;  and  k is  in  m^ . 

Index  Y in  Eq(31)  is  an  expression  for  the  ratio  of  the  expected  average  indoor 
radon  level  to  the  EPA  threshold  level  of  4 pCi/2.  Thus  if  Y > 1 it  is  expected 
that,  unless  special  precautions  are  taken,  the  EPA  threshold  level  will  be 
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exceeded.  Y is  calculated  from  the  invariant  index  properties  previously 
discussed,  namely:  the  dry  in-place  density  of  the  soil,  the  radium  activity 
concentration,  the  porosity,  and  the  dry  gas  permeability. 


Except  for  Eq(28)  which  considers  air  change  periods,  the  previous  expressions 

do  not  consider  the  characteristics  of  the  buildings,  even  though  such 

consideration  is  implicit  in  the  expressions  developed.  Other  equations  for 

source  potential,  however,  also  consider  parameters  associated  with  the 

building.  A calculated  potential  in  terms  of  activity  per  unit  time  (pCi-h"^  or 

Bq-s'^)  was  presented  by  Nazaroff  and  Sextro,  1989.  Source  potential  is  defined 

as  the  "maximum  sustainable  entry  rate"  based  on  a specified  indoor-outdoor 

pressure  difference,  isotropic  soil  characteristics,  dominance  of  convective 

transport,  and  air  penetration  at  the  perimeter  basement  wall  through  cracks  of 
an  assumed  widths: 


For:  p^k  < 0 . 5n/iH^  • ln(2H/r) 


F 




^ ln(H/r) 


...(32a) 


for  larger  values  of  k (approximately  k > 3. 5- 10”^^  m^)  : 

2/3 

...(32b) 


f H ] 

1 4A  J 1 

[ ln(H/r)  J 

where : 


H = depth  of  a perimeter  crack  between  the  basement  slab  and  the  wall 
below  the  grade  outside  the  house 

r = half  the  width  of  a perimeter  crack 

L = length  of  the  perimeter  crack. 


For  a house  with  a 40  m perimeter  crack  in  the  basement,  a 450  m^  volume  and  an 
air  exchange  rate  of  2 hours,  F should  not  exceed  10  Bq-s'^. 


Equations  32  a and  b were  developed  for  steady-state  flow,  which  account  for  the 
fact  that  the  radon  entry  rate  is  directly  proportional  to  k.  When  dimension  H 
is  interpreted,  account  must  be  taken  of  the  fact  that  the  air  flow  is  below  the 
wall  foundation.  It  was  suggested  in  Eq(19)  to  add  1.5m  to  H in  order  to  account 
for  this  fact.  A more  reasonable  approach  may  be  to  measure  H from  the  bottom  of 
the  footing. 
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Equations  32  were  purposely  written  in  a way  which  separates  the  terms 
associated  with  site  characteristics  from  those  associated  with  building 
characteristics.  These  terms  can  be  de-coupled  for  the  purpose  of  site 
characterization. 


Another  approach  to  source  potential  determination  was  developed  by  Tanner, 
1987.  A "mean  migration  distance"  is  calculated  by  the  following  equation: 


M = 


V 

2(nAD)®-5 


4nAD 


[(nD)A-l] 


..(33) 


V is  taken  such  that  p decreases  linearly  by  5 Pa  over  a distance  of  4.6-M/n 
which  is  the  distance  over  which  the  radon  concentration  would  decay  to  1%  if 
the  flow  velocity  is  as  in  Eq(33).  Solutions  are  obtained  by  iteration.  A 
volume  of  soil  is  calculated  by  folding  the  foundation  walls  outward  from  the 
building  and  multiplying  the  projected  area  by  M.  The  radon  contained  within 
the  pore  space  of  this  volume  is  asstomed  to  be  the  maximum  available  radon.  The 
difference  between  this  approach  and  approaches  previously  discussed  is  that 
diffusive  as  well  as  convective  flow  is  considered.  A potential  problem 
envisioned  is  that  the  soil  volume  is  a linear  function  of  k,  while  radon 
concentrations  observed  appear  to  be  approximately  proportional  to  k°-^.  The 
reason  for  this  dimensional  relationship  is  that  a steady-state  flow  condition 
(linear  pressure  gradients)  is  used  to  derive  M. 


As  in  Eq(31),  equations  32  and  33  can  be  expressed  in  terms  of  measurable 
invariant  soil  index  properties.  Assumptions  would  then  have  to  be  made  to 
account  for  the  effect  of  soil  moisture  content  on  permeability,  diffusion,  and 
radon  concentration  in  the  soil  gas. 

At  the  present  time  it  would  be  premature  to  recommend  one  of  the  expressions 
discussed  in  this  section  and  exclude  the  other  expressions  from  consideration. 
This  choice  should  be  made  after  comparison  with  actual  field  data  and  some 
further  parametric  studies.  Eqs.(28),  (29),  and  (31)  use  invariant  soil 
properties.  Eq  (30)  would  be  similar  to  Eq  (31)  when  modified  to  be  used  with 
invariant  soil  properties.  Eqs . (32)  and  (33)  could  also  modified  to  be  used 
with  invariant  soil  properties  with  an  allowance  for  soil  moisture  content. 
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Comparison  with  New-York  State  data  [Kunz  et  al.,  1988]  indicates  that  the  trend 
of  the  data  is  very  well  predicted  by  Eq  (31),  which  was  calibrated  to  fit  the 
average  value  of  the  data  (there  is  a good  fit  with  individual  data  points) . 

A source  potential  measure  based  on  invariant  soil  properties  will  have  to  take 
several  conditions  into  consideration,  probably  in  the  form  of  modifiers.  These 
include;  climatological  conditions  such  as  rainfall,  wind  intensities,  length  of 
winter  season  and  depth  of  frost  cover;  and  soil  profile  characteristics  such 
as  drainage  characteristics,  depth  to  rock,  soil  stratification,  location  of 
phreatic  surface  and  overconsolidation  of  the  deposits. 

5.  SITE  EXPLORATION  FOR  RADON  SOURCE  POTENTIAL 

5.1  Criteria 

Three  criteria  govern  the  procedures  considered  in  this  report, 

(1)  The  site  is  characterized  by  soil  index  properties  which  do  not  vary 
as  a result  of  transient  conditions. 

(2)  The  measurements  should  not  depend  on  subjective  judgment  (though 
their  interpretation  involves  professional  judgment). 

(3)  While  approximate  procedures,  such  as  manual-visual  soil 
identification  [ASTM  D2488]  or  estimates  of  permeability  by  Table  5 
are  desirable  for  practical  and  economic  reasons,  accurate  and 
repeatable  methods  such  as  ASTM  D2487  and  appropriate  insitu  and 
laboratory  permeability  tests  that  can  be  used  as  a backup  in  case  of 
doubt  or  disputes  must  be  identified. 

5.2  Options  for  Site  Exploration 

Several  options  are  outlined  here  on  a preliminary  basis.  These  options  are  not 
mutually  exclusive.  Neither  are  they  intended  as  a substitute  for  present  site 
exploration  methods  which  measure  in  situ  conditions  existing  at  any  given  point 
in  time.  Rather,  the  proposed  site  exploration  methods  are  intended  to  be  used 
to  provide  information  prior  to  construction  and  information  which  is 
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independent  of  site  conditions  which  may  be  temporary.  One  of  the  proposed 
methods  can  be  readily  integrated  with  the  Standard  Penetration  Test,  which  is 
presently  the  most  commonly  used  soil  exploration  procedures. 

5.2.1  Test  Borings 

The  Standard  Penetration  Test  [ASTM  D 1586]  is  one  of  the  most  common  soil 
exploration  methods  used  in  the  U.S.  Because  it  is  so  commonly  used,  it  will  be 
in  many  instances  practical  to  combine  its  use  for  a routine  geotechnical 
investigation  with  a determination  of  radon  source  potential.  The  test  is  also 
relatively  inexpensive  and  equipment  and  trained  personnel  are  readily  available 
in  any  U.S.  location. 

The  test  drives  a standardized  sampler  into  the  ground  with  a hammer  which 
delivers  a standardized  amount  of  potential  energy  with  each  blow.  A partially 
disturbed  sample  is  retrieved  which  has  a 1 3/8  in. (35  mm)  diameter  and  is  up  to 
18  inches  (457  mm)  long.  It  is  realistic  to  expect  that  a 5-inch  (127  mm)  long 
sample  can  in  most  instances  be  retrieved  without  difficulty.  This  provides  a 
sample  of  approximately  200  g dry  weight,  large  enough  for  radium  activity 
concentration  measurement  by  7 - spectroscopy.  The  sample  is  also  large  enough 
for  soil  classification  in  accordance  with  ASTM  D2487.  Normal  procedures 
include  a field  manual-visual  soil  classification  by  ASTM  D2488,  and  a 
laboratory  classification  by  ASTM  D 2487.  The  test  does  not  provide  an  accurate 
measurement  of  in  situ  dry  density,  because  the  sampler  is  driven  into  the 
ground,  partially  disturbing  the  soil.  Density  could  be  estimated  from  the 
driving  resistance.  Alternately,  the  height  of  the  sample  retrieved  from  the 
split-spoon  sampler  can  be  recorded  in  the  field.  Dry  density  can  then  be 
approximately  determined  in  the  laboratory.  After  completion  of  the  laboratory 
classification,  dry  samples  will  be  sealed  and  stored  for  30  days  for 
determination  of  radium  activity  concentration. 

The  test  will  provide  a depth  profile  as  required,  from  which  the  following 
information  can  be  extracted:  depth  to  rock  (if  needed),  location  of  phreatic 
surface,  n,  Aj^^ , k (estimated  from  particle  size  distribution  and  density), 
natural  moisture  content,  and  Atterberg  limits  (if  appropriate)  which  in  turn 
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can  be  used  to  determine  swelling  potential  (propensity  for  fissure  formation, 
refer  to  Table  4)  . It  is  also  suggested,  in  addition  to  the  Atterberg  limit 
test,  to  perform  a drying  test  [Yokel  and  Tucker,  1980]  in  order  to  supplement 
the  information  from  the  Atterberg  limits.  This  information  should  be  adequate 
for  an  approximate  determination  of  the  radon  index.  If  desired,  the  saturation 
(m)  can  be  estimated  from  the  natural  moisture  content. 

If  the  soil  is  finegrained,  more  accurate  data  can  be  obtained  by  retrieving  an 
essentially  undisturbed  thin-wall  tube  sample  in  accordance  with  ASTM  D1587. 
Such  a sample  will  permit  accurate  laboratory  determination  of  7^,  n and  k. 
Examination  of  the  sample  in  the  laboratory  can  also  provide  clues  on  the 
presence  of  fissures. 

More  recently  a "continuous  sample  tube  system"  has  been  introduced,  in  which  a 
3 inch  (76  mm)  O.D.,  2 1/2  inch  (64  mm)  I.D,  core  barrel  is  inserted  in  a 
hollow-stem  auger  to  extract  a 5 ft  (1.52  m)  long  relatively  undisturbed  sample. 
This  latter  method  permits  extraction  of  relatively  undisturbed  samples  of 
granular  soils  or  desiccated  clays.  The  sample  allows  a much  more  accurate 
determination  of  dry  density  than  the  SPT.  These  samples  may  also  be  used  to 
detect  fissures  in  fine  grained  soils  and  in  some  instances  for  laboratory 
permeability  testing. 

5.2.2  Manual  Procedure 

A simple  manual  procedure  is  also  considered,  particularly  when  an  isolated  site 
is  explored  where  it  would  be  expensive  to  bring  in  SPT  equipment.  The  suggested 
procedure  includes:  measurement  of  the  radium  activity  concentration  with  a 
portable  7-ray  spectrometer;  measurement  of  soil  density  and  moisture  content 
with  a portable  nuclear  gage;  manual-visual  soil  classification  in  the  field; 
and  retrieval  of  a soil  sample  (approximately  1 kg)  for  laboratory  soil 
classification.  This  procedure  would  provide  a rapid  and  relatively  inexpensive 
evaluation  of  radon  source  potential.  However,  readily  available  instruments 
would  be  mainly  usable  for  surficial  measurements  (unless  a test  pit  is  dug) . It 
is  proposed  to  build  portable  down-hole  spectrometers  and  density  gages,  which 
can  be  used  in  an  auger  hole.  Such  instruments  have  been  developed,  but  they  are 
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not  readily  conraiercially  available  and  their  portability  may  have  to  be 
improved. 

5.2.3  Percolation  Test 

Percolation  tests  [HUD  Minimum  Property  Standards,  1973]  which  are  routinely 
required  for  many  rural  building  sites  could  be  utilized  to  estimate 
permeability  and  integrated  with  the  manual  procedure  in  2.  This  would  improve 
the  accuracy  of  the  estimate  of  k.  However,  a study  is  required  to  correlate 
test  results  with  k. 

5.2.4  Existing  Radon  Testing  Protocols 

Existing  radon  testing  protocols  [i.e.  Tanner,  1987]  can  be  used  to  measure 
radon  concentration  in  the  soil  gas  and  gas  permeability,  together  with 
determinations  of  in  situ  moisture  content  and  dry  density.  The  results  can 
than  be  corrected  for  moisture  content  to  obtain  an  estimate  of  invariant  index 
properties.  Cm  can  be  measured  at  a great  enough  depth  to  provide  a measure  of 
r 

max  ' 

5.2.5  Backup  procedures 

In  some  instances,  when  the  accuracy  of  field  test  interpretations  is  doubted  or 
disputed,  a more  accurate  procedures  may  be  required  to  determine  k.  Several 
backup  procedures  are  available.  If  the  soil  is  dry  or  moist  a soil  gas 
permeability  test  (i.e.  Tanner,  1987)  can  be  conducted  with  an  estimated 
correction  factor  for  saturation.  In  wet  soils  this  procedure  would  be 
inaccurate.  Alternately,  k can  be  determined  by  ASTM  D2434  (constant  head  test) 
ASTM  D2435,  one-dimensional  consolidation  test,  or  horizontal  capillary  test 
[Lambe,  1958].  Water  permeabilities  will  have  to  be  corrected  for  the 
Klinkenberg  effect  (the  correction  factor  ranges  approximately  from  1 for  sands 
to  2 for  silts) . The  horizontal  capillary  test  results  must  additionally  be 
corrected  for  partial  saturation  (a  factor  ranging  from  2 to  3). 
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5.3  Selection  of  Exploration  Procedure 


The  exploration  procedure  selected  in  any  particular  instance  would  depend  on 
several  conditions.  If  a single,  isolated  housing  site  is  involved,  the  manual 
procedure  is  probably  the  most  practical,  because  it  is  expensive  to  mobilize 
the  drilling  equipment.  The  manual  procedure  would  also  be  economical  when 
percolation  tests  have  to  be  conducted,  because  the  two  procedures  can  be 
combined.  On  the  other  hand,  if  many  housing  sites  are  involved  there  would  be 
an  economy  of  scale  and  the  test  boring  procedure  would  be  more  practical.  The 
test  boring  procedure  is  also  an  attractive  and  inexpensive  option  when  the 
radon  source  potential  test  can  be  combined  with  a geotechnical  exploration 
required  for  other  purposes. 

The  test  boring  procedure  has  the  advantage  of  providing  a depth  profile.  While 
the  manual  procedure  can  also  provide  a depth  profile,  it  may  become  expensive 
when  the  desired  depth  of  exploration  exceeds  1.5  m. 

6 . STANDARD  DEVELOPMENT 

This  study  indicates  that  it  is  feasible  and  practical  to  develop  testing 
protocols  for  the  determination  of  invariant  index  properties  of  the  subsoil 

associated  with  radon  source  potential.  However,  it  is  necessary  to  establish 
the  credibility  of  a standard  for  the  determination  of  the  radon  source 
potential  of  sites  which  is  based  on  these  invariant  soil  properties.  The 

following  program  is  proposed  to  accomplish  this  goal. 

Phase  I - Feasibility  study:  This  stage  is  completed  with  this  report. 

Phase  II-  (a)  Development  of  testing  protocols 

(b)  Field  and  laboratory  study  (effect  of  moisture  content  on  gas 
permeability,  comparison  of  water  and  gas  permeability,  study  of 
HUD  percolation  test,  development  and  tryout  of  instrument 
packages) 

(c)  Parametric  study  (use  transient  model  to  assess  effect  of 

various  parameters  and  examine  existing  source  potential 

equations . 
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Phase  III-  Data  acquisition  and  verification  of  predictive  models 
Phase  IV  - Preparation  of  draft  standard 

Development  of  Testing  Protocols: 

An  outline  of  options  is  provided  in  Section  5.2.  These  options  must  be  spelled 
out  in  more  specific  terms  and  procedures  must  be  specified  to  avoid  ambiguities 
and  insure  the  quality  control  of  the  planned  data  acquisition. 

Field  and  Laboratory  Studies: 

Effect  of  moisture  content:  Present  data  on  the  effect  of  saturation  on  air 
permeability  are  inadequate.  A well  controlled  laboratory  study  is  needed  to 
provide  a data  base  similar  to  that  which  exists  for  moisture  effects  on  the 
diffusion  coefficient. 

Study  of  Klinkenberg  effect:  dry  gas  and  saturated  water  permeabilities  need  to 
be  compared  for  different  soil  types  to  provide  a reliable  data  base  for 
converting  water  permeabilities  to  gas  permeabilities  at  atmospheric  pressure. 

HUD  Percolation  Test:  No  study  could  be  located  to  correlate  the  test  with 
permeability.  Also,  specifications  for  the  test  should  be  tightened.  A similar 
test  was  extensively  used  by  the  Bureau  of  Reclamation  [Bureau  of  Reclamation, 
1963]  and  can  be  correlated  with  permeability,  however  the  latter  test  uses  a 
cased  hole  and  a larger  water  head  and  requires  a drill  rig.  It  is  not  clear 
whether  the  cost  of  the  test  would  substantially  increase  if  the  BUREC  test  were 
used.  A study,  including  field  data  acquisition,  and  an  appropriate  analytical 
model  for  the  test  is  needed. 

Development  of  Instrument  Package:  a portable  downhole  spectrometer  and  a 
downhole  nuclear  moisture-density  gage  are  needed.  These  instruments  would  not 
only  be  useful  for  the  manual  procedure  outlined  in  5.2,  but  they  could  also  be 
used  together  with  the  Tanner  protocol  and  in  conjunction  with  the  proposed  SPT 
procedure . 
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Parametric  Study 


It  is  proposed  to  use  a computer  code  combining  transient  convective  flow  with 
diffusion  and  radon  generation  to  test  the  source  potential  models  described  in 
Section  4 and  if  appropriate  develop  a new  model.  The  calculation  will  converge 
to  steady- state  convective  flow  if  the  duration  of  the  pressure  drop  is  long 
enough . 

Data  Acquisition  and  Verification 

The  testing  protocols  developed  in  Phase  II  will  be  used  to  determine  the 
relevant  invariant  soil  parameters  for  sites  where  extensive  radon  data  are 
available.  On  some  of  the  sites,  presently-used  in  situ  exploration  protocols 
(i.e.  Tanner,  1987)  will  be  used  along  with  the  testing  protocols  developed  in 
Phase  II  to  obtain  comparative  data.  The  sites  selected  will  include  sites 
where  radon  measurements  are  presently  performed  by  Federal  agencies  and  other 
sites  where  data  are  available.  An  initial  quick  pilot  study  will  be  conducted 
to  insure  that  the  data  sample  will  cover  a sufficiently  broad  range  of  soil 
permeabilities  and  other  characteristics. 

The  results  of  the  study  will  be  plotted  against  various  equations  for  source 
potential  to  determine  how  well  they  predict  the  source  potential  in  the  data 
sample.  A best  fit  source  potential  equation  will  be  developed. 

Preparation  of  Draft  Standard 

A draft  standard  for  site  characterization  will  be  proposed.  The  standard  will 
include  alternate  testing  protocols  and  methods  of  determining  source  potential. 
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